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Abstract 

We use the integrable Kaup-Boussinesq shallow water system, modified by a small 
viscous term, to model the formation of an undular bore with a steady profile. The 
description is made in terms of the corresponding integrable Whitham system, also ap- 
propriately modified by friction. This is derived in Riemann variables using a modified 
finite-gap integration technique for the AKNS scheme. The Whitham system is then 
reduced to a simple first-order differential equation which is integrated numerically to 
obtain an asymptotic profile of the undular bore, with the local oscillatory structure 
described by the periodic solution of the unperturbed Kaup-Boussinesq system. This 
solution of the Whitham equations is shown to be consistent with certain jump condi- 
tions following directly from conservation laws for the original system. A comparison 
is made with the recently studied dissipationless case for the same system, where the 
undular bore is unsteady. 



Undular bores are nonlinear wave-like structures, which are generated in 
the breaking profiles of large-scale nonlinear waves propagating in dispersve 
media. A general theory, based on the Whitham modulation equations, has 
been previously developed for dissipationless, unsteady, undular bores on the 
basis of completely integrable models such as the Korteweg-de Vries equation, 
nonlinear Schrodinger equation etc. The introduction of physically important 
small dissipation in the system dramatically changes its properties, allowing 
in some cases for the presence of steady solutions. The most explored model 
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describing the effects of friction on an undular bore is based on the uni-directional 
Korteweg-de Vries equation, modified by a small friction term, which can take 
various forms. Appropriate perturbation techniques have been used to obtain 
asymptotic solutions. However, unlike the case for conservative undular bores, 
no general approach seems to be available. Here, using an integrable version 
of the bi-directional Boussinesq equations, but modified by a small Burgers- 
like dissipation term, we develop a modulation theory of frictional shallow water 
undular bores, which can also be extended to other non-conservatively perturbed 
integrable systems. 

1 Introduction 

It is well known that solution to an initial value problem for the inviscid dispersionless shallow 
water equations may develop a wave-breaking singularity after a finite time, when the first 
spatial derivatives become infinite. After the wave-breaking point, a formal solution becomes 
multi- valued and loses its physical meaning. The divergence of the spatial derivatives at the 
wave-breaking point suggests that higher-order terms must be taken into account. These 
terms can be either dispersive or dissipativc in nature, or, as here, a combination of both. 
The form of the solution after the breaking time then strongly depends on the comparative 
values of the dispersion and dissipative terms. If dissipation can be neglected in favour 
of dispersion, the solution in a certain neighbourhood of the breaking point assumes the 
form of an expanding nonlinear oscillatory structure with a solitary wave train generated in 
the vicinity of its leading edge. This structure provides a dispersive resolution of a breaking 
singularity, and is an unsteady undular bore (or a dispersive shock in a different terminology). 
Unsteady undular bores have been studied extensively in the last thirty years on the basis 
of exactly integrable nonlinear wave equations. The original formulation of the problem was 
given by Gurevich and Pitaevskii (1973) who proposed to describe the expanding coUisionless 
shocks (a plasma analog of undular bores) with the aid of the Whitham-avcragcd equations 
for the integrable Korteweg-de Vries (KdV) equation. The Gurevich-Pitacvskii theory has 
been extended to other important integrable systems such as the nonlinear Schrodingcr and 
Kaup-Boussinesq equations (see Kamchatnov (2000), for instance, for the detailed account 
on the Gurevich- Pitaevskii theory). 

The introduction of small dissipation can, in some cases, balance the dispersive effects 
so that the undular bore eventually acquires a steady profile, but remains oscillatory in 
space. The analytical study of steady (frictional) undular bores was initiated in the classical 
work of Benjamin and Lighthill (1954) on shallow water waves. Another important work on 
the same subject, but in the context of coUisionless plasma shocks with small dissipation, 
is Sagdeev (1964). In both works the authors use a mechanical analogy with a weakly 
damped nonlinear oscillator to explain the main observable features of undular bores: the 
formation of the lead solitary wave and degeneration into linear sinusoidal waves at the rear. 
It was also suggested that the undular bore transition conditions must be consistent with 
the conservation of mass and momentum across the transition zone, while the violation of 
hydrodynamic energy conservation is remedied by taking into account the generated waves. 

A simple model with an analytic description of shallow water frictional undular bore 
is provided by travelling wave solutions of the KdV-Burgers equation (see for instance 



2 



Whitham(1974)) 

Ut + UUr, + U^xx = l^Uxx (1) 

with a small dissipation coefficient < u <^ 1. A detailed study of such solutions was 
made by Johnson (1970), who applied a direct perturbation procedure (Kuzmak 1959) to 
the periodic solution of the unperturbed {u = 0) equation (the KdV cnoidal wave), and 
performed matching of the leading order approximation with the solitary wave to obtain 
a closed description. Johnson's (1970) solution has been used by Smyth (1988) for the 
description of the effect of small dissipation on resonant flow over topography. 

The description of the undular bore on the basis of the steady travelling wave solutions of 
the KdV-Burgers equation has two inherent restrictions: (i) as it is based on a unidirectional 
equation, it does not reveal the transition (jump) conditions across the undular bore and 
can link any two given constant states u = U2 and u = Ui, U2 > Ui, (ii) it describes only the 
established (steady) regime, and says nothing about the undular bore formation. The further 
development of the Whitham modulation theory (Whitham (1965, 1974)) in 1970-80s due to 
Gurevich and Pitaevskii (1973, 1987), Lax, Levermore and Venakides (see the review (1994) 
and references therein), Flaschka, Forest and McLaughlin (1982), Dubrovin and Novikov (see 
the review (1989) and references therein) and many other authors, has made it clear that the 
consistent description of undular bores (both conservative and frictional) should be made 
in the framework of the hydrodynamic-type Whitham equations describing the evolution of 
nonlinear modulated waves. Although the Whitham equations, based on averaging over the 
periodic wave family, correspond to the leading order of a direct perturbation procedure, 
which formally diverges when the wave period tends to infinity ( the solitary wave limit), 
their solutions reveal only a weak singularity at the leading edge of the undular bore (see 
Gurevich and Pitaevskii (1973, 1987)) and yield the correct value for the lead solitary wave 
amplitude (while its position, of course, is not determined accurately). Also, the Whitham 
equations have been shown to inherit an integrable (or perturbed integrable) structure from 
the original system and allow in some cases the effective construction of exact global solutions 
using powerful methods developed in the theory of finite-gap integration, and in the theory 
of integrable Hamiltonian systems of hydrodynamic type (Tsarev (1985, 1990), Dubrovin 
and Novikov (1989)). 

The modulation theory of the "integrable" shallow water undular bore has been con- 
structed by El, Grimshaw and Pavlov (2001) using the extension to the case of the bi- 
directional Kaup-Boussinesq system of the original formulation of Gurevich and Pitaevskii 
(1974) for the decay of an initial discontinuity in the KdV equation. A more general case of 
the formation of an undular bore in the vicinity of a "cubic" breaking point has been studied 
in (El, Grimshaw and Kamchatnov, 2005). An asymptotic theory of the formation of soliton 
trains from a "big" enough initial pulse for the Kaup-Boussinesq system was developed by 
Kamchatnov, Kraenkel and Umarov, 2003. 

The modulation equations for the KdV-Burgers equation were derived by Gurevich and 
Pitaevskii (1987) and Avilov, Krichever and Novikov (1987). Other forms of the dissipa- 
tive term have been considered by Gurevich and Pitaevskii (1991) (boundary layer-type), 
and Myint and Grimshaw (1995) (boundary layer dissipation and Rayleigh friction). In 
their work on the modulation theory of the KdV-Burgers equation, Gurevich and Pitaevskii 
(1987) deduced an exact steady solution of the Whitham system corresponding to the steady 
undular bore. It is worth noting that their solution exactly coincides with the leading-order 
perturbation solution to obtained much earlier by Johnson (1970). Avilov, Krichever and 
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Novikov (1987) have shown numerically that this solution is indeed the large-time asymptotic 
modulation solution of (Q) with the initial conditions in the form of a smooth step. 

A general procedure for obtaining perturbed modulation system for the KdV equation, 
based on the finite-gap integration machinery, was formulated by Forest and McLaughlin 
(1984). A more effective method for the case of periodic modulated waves, applicable not 
only to the perturbed KdV equation, but also to the whole perturbed AKNS hierarchy has 
recently been designed by Kamchatnov (2004). 

In this paper, we apply Kamchatnov's procedure to the Kaup-Boussinesq system mod- 
ified by a small Burgers-like dissipative term. Being bi-directional, this system allows a 
more realistic modelling of frictional undular bores than the KdV-Burgers equation (^. We 
distinguish several characteristic stages of the undular bore evolution and construct exact so- 
lutions of the Whitham equations describing the initial (unsteady) and final (steady) stages 
of evolution of the KBB undular bore. The methods used in this paper also allow for the 
analytic description of undular bores generated in presence of an external force. In partic- 
ular, an important class of problems of this kind occurs in the description of resonant flow 
over topography (see Grimshaw &; Smyth (1986) and Smyth (1987, 1988)). 



2 Formation of a frictional undular bore: general de- 
scription 

We consider formation of a frictional undular bore in the Kaup-Boussinesq system modified 
by a small viscous term. In dimensionless variables this system has the form: 

ht + {hu)x + \u^^^ = 0, .^s^ 

Ut + UUx + hx = fUxx, 

where h{x,t) denotes the height of the water surface above a horizontal bottom, u{x,t) is 
related to the horizontal velocity field averaged over depth (see (Kaup (1976) for the detailed 
derivation of the inviscid system) and < z/ ^ 1 is a small viscosity coefficient. We shall 
call Eq. Q the Kaup-Boussinesq-Burgers (KBB) system. 

Note that the frictional term appears only in the second equation, which represents the 
momentum balance, and is absent in the first equation which represents the mass balance. 
Also, the derivation of this system requires that the frictional term is a small term, of the 
same order as the small dispersion term. In the sequel, however, we will be treating the 
frictional term as a small perturbation to an inviscid system. 

Compared to the KdV-Burgers equation (^, the KBB system has the essential ad- 
vantage of modelling bi-directional wave propagation, so the undular bore description would 
necessarily include transition conditions, which should be consistent with the jump condi- 
tions following from the conservation laws of the system On the other hand, the KBB 
system (j2)) is a perturbed integrable system, which retains the advantage of amenability to 
an effective analytic study. A drawback of the KBB system as a model system is the presence 
of a high-wavenumber instability of the constant solutions. That is, the linearized equations 
allow for growing waves at large wavenumbers (see El, Grimshaw, & Pavlov 2001). 

Also, there is the disadvantage for the description of undular bores that there is no 
"physical" momentum conservation law, which leads to formally "nonphysical" transition 
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conditions. We will show, however, that the transition conditions following from the solu- 
tions of Eq. (j2)) are asymptotically consistent with the classical jump conditions for shallow 
water bores within the range of applicability of the KBB system. 

We consider initial data at t = to: hQ{x) = h{x,tQ), uq{x) = u{x,to) for the system Q in 
the form of a smooth transition between two constant states: 

h = hi, u = Ui as x ^ +00, , . 

h = h2, u = U2 as X ^ —00, 

so that /12 > hi and the characteristic width of the transition region / ^ 1. There 
are two typical spatio-temporal scales associated with the initial- value problem 0, (jSI): 
characteristic "nonlinear-dispersive" scale At„;j„ ~ Axnim ~ 1 and the "dissipative" scale 
Atd ~ Axd ~ z/^^. Let us suppose that the initial conditions are chosen such that td 3> tnZm, 
i.e. we have 

/ > 1, < 1, and z// < 1 . (4) 

Then, following Avilov, Krichever and Novikov (1987), we distinguish several stages in 
the process of the formation of a frictional undular bore, and discuss some limitations of the 
applicability of this scenario. 

Stage 1. to < t < t^r ^ where the breaking time tbr will be defined below; ttr — to ~ I- 
Due to Eq. Q the initial data satisfy the following inequalities: 

^\uo"'\ < lihoUoYl, z/|uo"| < l^o^iol, (5) 
Thus this stage of the evolution can be described by the ideal shallow water system 

ht + {hu)n, = 0, ut + uun, + /la; = 0, (6) 
which can be represented in the diagonal form 

d\+ 1 , , , , 9A+ (9A_ 1 , , , , dX- , . 

Here 

A± = I ± (8) 

are the Riemann invariants of Eqs. (jH)). 

Initial data are given by two functions A-t-(x, to) and A_(x,to) determined by the initial 
distributions hQ{x) and uq{x). The system has two families of characteristics in the (x, t) 
plane along which one of two Riemann invariants (either A4. or A_) is constant. The wave- 
breaking point corresponds to the moment when characteristics of one of the families begin 
to intersect, so that the corresponding Riemann invariant becomes a three- valued function 
in the physical plane. Let such an intersection occur for the characteristics transferring 
the values of A_|_. Then at the wave-breaking point the profile of A+ as a function of x 
has a vertical tangent line and, hence, in vicinity of this point it varies very fast, whereas 
the second Riemann invariant varies with x much slower and may be considered here as a 
constant parameter: 

A_ = Ao = const. (9) 
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Thus, in the vicinity of the breaking point at t = t^r we are deahng with a simple wave. 
Without loss of generality we choose tf,r = 0. The second equation in ((Tj) is identically 
satisfied by Eq. Q. The first equation in ((7j) then has the well-known solution 

x-i(3A+ + Ao)t = /(A+), (10) 

where /(A+) is an inverse function to an initial profile A+(x, 0). Due to our normalization, 
the function x = /(A_|_) must have an inflexion point with a vertical tangent line at t = 0. 
In the vicinity of this point /(A+) can be approximated by a cubic function, 

X - ^(3A+ + Xo)t = -C(A+ - A° )3, (11) 

where C and A^ are constants. Without loss of generality Eq. ()1H) can be cast into the form 
(see El, Grimshaw and Kamchatnov (2005) for details) 

x-^(3A+ + Ao)t = -A^. (12) 
It corresponds to the wave breaking picture shown in Fig. ^ 




Figure 1: Wave breaking of the water elevation in the dispersionless limit; Aq is taken equal 
to -10. 

Stage 2. Ur < t <^ v'^ . At this stage, dispersion should be taken into account in the 
vicinity of the breaking point which implies consideration of the integrable KB system: 

ht + {hu)^ + ^Uxxx = 0, ut + uux + = 0. (13) 

with the initial data following from (©, ©: 

t = 0: ^ + v^=-xi/^ ^-v^=Ao. (14) 

The combined action of nonlinearity and dispersion leads to the generation of an expanding 
nonlinear oscillatory structure occupying the finite space interval (x~; x^) (see Fig. |2)). This 
structure is an unsteady, "conservative" undular bore. Outside the interval {x~;x~^) the 
flow is smooth and is described by the solution ©, (fT^. The solution of the problem now 
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consists of two parts. Following Gurevich and Pitaevskii (1973), we suppose that the region 
of oscillations can be approximated by a modulated periodic solution of the KB system. 
Its global evolution is then determined by the Whitham equations and the problem reduces 
to finding the solution of the Whitham equations that matches the solution ()12|) at the 
end points of the oscillatory region. One may say that this oscillatory region (the expanding 
undular bore) "replaces" a non-physical multi- valued region of the solution (fT^ . One should 
emphasize, however, that the boundaries of the undular bore do not coincide with the 
boundaries of the formal multi- valued solution. 
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Figure 2: Initial stage of the undular bore development. The plot corresponds to the time 
t = 1 and Ai = —10. Dashed line shows the formal solution in the dispersionless limit. 



The corresponding modulated solution of the KB system has been constructed by El, 
Grimshaw and Kamchatnov (2005). Here we briefly outline the resulting formulas. The 
derivation of the complete modulation system with the account of the dissipative corrections 
will be presented in Section 4. 

The local wave form of the undular bore is given by the periodic travelling wave solution 
of the KB system (fT^ . and is given by the expressions 

m(x, i) = Si- 2^{6), h{x, t) = \sl - S2 - 2/^^(6*) + Si/i(6'), 6 = x- |sit, (15) 

where 

A2(A3 - Ai) - Ai(A3 - A2)sn2 f V(A4 - A2)(A3 - Ai) ^, m) 

m = } = (16) 

A3 - Ai - (A3 - A2)sn2 (V(A4-A2)(A3-Ai)^,m] 

A4 ^ A3 ^ A2 ^ Ai. 
Here sn{6, m) is the Jacobi elliptic function with the modulus 

(A4-A,)(A3-A,)- * ' 

The connection of the constants si, S2 in Eq. p5|) with the parameters \j in Eq. (fT^ is given 
by 

4 

Si = ^ Aj , ■52 = ^ . (18) 
j=l i<j 



The soliton limit (m = 1) is obtained either for Ai = A2 or for A3 = A4. For A3 = A4, which 
corresponds to right-propagating sohtons, Eq. (fTH|) yields 

A2 - Ai + (A4 - A2)/cosh2[v/(A4 - A2)(A4 - Ai) 6] 

Modulations Xi{x, t) in the travelling wave solution are described by the Whitham equations, 
which have been derived for the KB system by El, Grimshaw, and Pavlov ( 2001) in the 
Riemann form (see also El, Grimshaw and Kamchatnov (2005)). 

^ + t;,|^ = 0, z = 1,2,3,4. (20) 
Here the characteristic velocities Vi are expressed in terms of Xj as 

^-^*)''- ^'"l-' (21) 

where the phase velocity V and the wavelength L are given correspondingly by 

V^l±X., ' r , (22) 

K(m) being the complete elliptic integral of the first kind. 

The solution of the Whitham system ()20|) matching the dispersionless solution Q, (fT^ 
is obtained using the generalized hodograph transform (Tsarev 1985, 1990) and has the form 



+ _ 16„,,(3) I 8 \ „,,(2) , 2 \ 2 (1) I 1 \ 3 „• _ o q /i. 

Ai = Ao = const. 



where 

wl"^ = (l - ^^«) ^^'^ ^ = 1' 2' 3, 4. (24) 
Here the functions W^''\Xi, . . . , A4) are obtained as coefficients of the series expansion 

W = I ^ ^ =J2^ = i + lsrl+ilsl - - ls,S2 + 1^3)-^+. . . . 

n(A-A,) 

(25) 

One can see that wl^'' = Vi coincide with the characteristic velocities ()21|). 

Formulas (f23|) - (j25|) define A2, A3, A4 implicitly as functions of x and t and, together with 
the travelling wave solution p5|l - (fT?)|l . determine the evolution of the undular bore at stage 
2. Dependence of the Riemann invariants on x at some fixed moment of time is shown in 
Fig.H 

The dynamics of the edges x^{t) of the undular bore at this stage 2 is given by the 
formulas 

x^{t) = \x,t+\^-ei\ (26) 




Figure 3: Dependence of Riemann invariants A2, A3, A4 on x at fixed time t = 1 and with 
Ai = —10. The dashed hne shows the corresponding dependence of A+ for the formal multi- 
valued solution of the KB equations in the dispersionless limit. 



x-(t)-iAot-^t^/^ + g^, V3t«|Ao|. (27) 

We note that this solution is generically realized only at the initial stage of the undular 
bore development, where the cubic approximation ()14|1 of the initial function is valid. After 
that, one should use the solution of the Whitham equations corresponding to the actual 
initial data. Such a solution can also be constructed in a closed form using the generalized 
hodograph method (see Gurevich, Krylov, El (1992) for the KdV case and El and Krylov 
(1995) for the defocusing NLS case). 

Undular bore 




Figure 4: Formation of an unsteady undular bore and rarefaction wave as an intermediate 
asymptotic for / < A(t) < z/~^ 

Stage 2a. If the dissipation coefficient u is small enough such that for some time interval 
the following inequality holds: / <^ ^{t) ^ ^^^j where A(t) = x^{t) — x~{t) is the undular 
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bore width, then an intermediate similarity asymptotic for the undular bore is reahzed where 
all the modulation parameters Xj in (|2Up are functions oir = x/t only. These solutions have 
been studied in detail in (El, Grimshaw and Pavlov 2001) where the problem of the decay of 
an initial discontinuity for the KB system has been considered. Generally, owing to the two- 
wave nature of the KB system such a solution would involve a rear rarefaction wave along 
with the leading undular bore (Fig. HJ. The similarity solution of the Whitham equations 
()20|1 in the undular bore region has the form: 

t;3 = f; Ai = aL') A4 = aJ\ (28) 

while the rarefaction wave is described by the similarity solution of the ideal shallow water 
equations ((Tj): 

1(A+ + 3A_) = |, A+ = aJ) (29) 

in the interval where 

aL'^ ^ A_ ^ aL'^ . (30) 

so that it matches a plateau region {A+ = A*^"*; A_ = A^"*} at its leading edge and the 
boundary values {A+ = X^^; A_ = A^^} at the trailing edge (see the diagram in Fig. El). 
Here 

Ai" = |±^. a£' = |±^, (31) 

It follows from the solution and relation (jSJ that two given constant states 

{hi.ui) and (/i2,U2), /i2 > hi could be connected with the aid of a single dissipationless 
undular bore (i.e. with no rarefaction wave generated), provided the following condition is 
satisfied 

y - V ^2 = Y - V ^1 • (32) 

One can notice that this condition (transition relation) coincides with the relationship be- 
tween flow parameters at any two points in the formal simple wave solution for the ideal 
shallow water equations ©. At the same time, as we stressed before, the solution of the 
Whitham equations does not coincide with the three-valued simple wave solution of the 
shallow water equations. 

We emphasize that the similarity stage of the undular bore evolution may not be realized 
at all if the dissipation coefficient is not small enough (see discussion in Avilov, Krichever 
and Novikov (1987) for the KdV-Burgers case). 

Stage 3. t z/"^. At this stage, the dissipation effects are accumulated to the degree 
that they begin to compete with the combined action of nonlinearity and dispersion. The 
dynamics of the undular bore is governed now by the full KBB system (0). The local wave 
form in the undular bore is still described by the periodic solution (jl5|) . (jl6p but the Whitham 
equations now become inhomogeneous 

dXi dXj , , 

= i= 1,2,3,4. (33) 

Explicit expressions for the functions Pi(Ai, . . . , A4) will be derived in Section 4 of this paper. 
We note that the undular bore at this stage is still unsteady. 
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Figure 5: Riemann invariant behaviour in the similarity asymptotic of the modulation solu- 
tion. 

Stage 4- t ^ v^^ . At this stage, the undular bore is reaching its steady regime so that it 
propagates as a whole with a single velocity, say c. The corresponding solution Aj = \j{x—ct) 
of the perturbed Whitham equations will be constructed in Section 5. Nevertheless, 
some important relationships for the steady frictional undular bore can be obtained from 
very general reasoning. 

Indeed, the original procedure for the derivation of the Whitham equations (Whitham 
1965) implies averaging the conservation laws of the original system over the periodic solu- 
tions. In the case of non-conservatively perturbed systems (which is the case for the KBB 
system) the averaging is performed over the periodic family of the unperturbed KB system 

In our case, we have only two conservation laws for the KBB system Q at our disposal: 

ht -f {hu + lu^x)x = 0, ut + {\u^ + h - vux)x = 0. (34) 

Of course, these are the KBB system Q itself. However, while the first equation in (33) is 
the representation of conservation of mass, the second equation is not that for momentum 
conservation, as this should have the term {hu)t and not Ut- Instead, the second equation is 
in effect the conservation of the Bernoulli expression. 

Averaging Eqs. over the periodic solution u{9), h{0) of the unperturbed KB system 
we obtain two modulation equations 

ht + (Hx = 0, ut+ (i^ + h)x = 0. (35) 

We note that the dissipative term drops out of the averaged conservation laws (jHKjl so the 
dissipation can only enter other modulation equations. But if the full perturbed modulation 
system admits the travelling solutions of the form f{x — ct), c being constant, then the 
equations must also admit such a solution. Substitution of /i = h{x — ct),u = u{x — ct) 
into Eqs. (jHHjl yields 

-ch + hu = -A, -cu + lv^ +h = B, (36) 
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A, B being constants. Let the established undular bore satisfy the boundary conditions Q 
at infinity. Then considering Eqs. (j36|) at x — ct — ±00 we obtain 

h2U2 - hiUi = c(/i2 - hi) , ^{ul - ul) + h2- hi = c{u2 - Ml) (37) 

which can be conveniently represented as 



c = ui + h2\ - — — -, U2 = ui + {h2- hi)^ - — — -. (38) 
y hi + h2 y hi + h2 

Thus, we have obtained an important restriction on the admissible family of the initial steps 
that may be eventually resolved into a single frictional undular bore with no additional 
(rarefaction) wave involved (cf. analogous condition ()32|) for dissipationless case). These 
conditions agree with the formal jump conditions obtained from the same two conservation 
laws for h and u of the ideal shallow water dynamics (jH} . 

However, it is well known that the usual "physical" jump conditions providing the 
mass and the momentum balance across the bore have the form (Benjamin, Lighthill 1954, 
Whitham 1974) 



c = u^ + h2^^^, «2 = «i + (^2-My%5^- (39) 

The discrepancy between the jump conditions (jHHj) and occurs due to absence of the 
proper momentum conservation law for the KB-Boussinesq system. This apparent disagree- 
ment, however, can be resolved by noticing that considered within the range of physical 
applicability of the KBB system, i.e. for small h2 — hi <^ hi the transition conditions (|HS)) 
and (j39p are asymptotically equivalent. In both cases we have 

r— 3 h2 — hi h2 — hi . ^ . 

C^Ml+V^l + 7 7i^5 M2 ~ Ml H 7=^- (40) 

4 y/hi y/hi 



In the next section we will show that the transition conditions in the form ()38|) also follow 
from the exact (non-periodic) travelling wave solution of the full KBB system (j21). 

3 Travelling wave solution of the KBB system: steady 
undular bore 

Here we shall study a steady travelling wave solution of the KBB system (j2)), i.e. we introduce 
the ansatz 

u = u{e), h = h{e), e = x-ct. (41) 

Its substitution into ^ leads to equations, which can be readily integrated once to give 

-ch + hu + \uee = -A, 

I 1 2 I T, I D ^ ' 

—cu + 2M + h = UUg + B, 

where A and B are again integration constants. Then the boundary conditions yield the 
relations and also 

^ _ hih2{u2 - Ml) 

h2-hi ' (43) 
B = Ui^Ui - c) + hi = U2{lu2 - c) + h2 
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U(u) 




u 



Figure 6: Potential ()45|) for an effective particle motion according to the Newton equation 




u 



Figure 7: Phase plane for an effective particle motion in the potential P5|l . 

Upon eliminating h from Eqs. one arrives at the equation 

uee + 4z/(m - c)ug = 4m(|m - c){u - c) - AB{u - c) - AA . (44) 

This can be considered as the Newton equation for a motion of a particle with a coordinate 
u (with 6 playing the role of time) in the potential 

U{u) = -\u^ + 2cM^ - 2{(? - B)u^ - A{Bc - A)u + constant (45) 

whose plot is shown in Fig.|Hl A phase plane (m, uq) for the undamped oscillator correspond- 
ing to the potential is shown in Fig. 7. There are three critical points (mi,0), (m2,0), 
(m3,0): the point (m2,0) is stable and the points (mi,0) and (m3,0) are unstable. The closed 
trajectories around the centre (m2, 0) correspond to a periodic motion and the separatrix 
corresponds to a soliton. Introducing small damping (the second term in the left-hand side 
of the Eq. (HD)) leads to an aperiodic oscillatory solution with the phase trajectory starting 
from the saddle point (mi,0) and eventually arriving after spiralling at the stable focus at 
(m2,0). This trajectory corresponds to a steady undular bore. 

The spatial oscillatory structure implied by this phase trajectory is the following. The 
large amplitude oscillations starting at (mi,0) correspond to the soliton train at the leading 
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edge of the undular bore and the small amplitude oscillations in the vicinity of (m2, 0) cor- 
respond to its trailing edge. It should also be noted that the configuration of the potential 
in Fig. 7 corresponds to the undular bore moving to the right. A bi-directional KBB system 
allows also an alternative configuration of the potential curve with the double roots at 
(ms, 0). This would reverse the picture so that the phase trajectory would start at the saddle 
point (^3,0) and after spiralling would again arrive at the potential minimum at (m2,0). 
This trajectory corresponds to the left-propagating undular bore 

The oscillatory profile of the bore can be found by numerical integration of Eq. 
However, a quite effective analytical theory can be developed on the basis of the Whitham 
modulation theory. The idea of the Whitham description of the frictional undular bore 
is to replace a weakly aperiodic motion of the particle in a given fixed potential with the 
asymptotically equivalent conservative motion in the potential which is slowly deformed. 
The Whitham equations then describe equivalent slow deformations of the potential. An 
advantage of the Whitham description in the case of perturbed integrable dynamics is that it 
utilizes the underlying integrable structure, and allows us to obtain the modulation equations 
using a universal technique based on powerful methods from finite-gap integration theory. 
At the same time, a straightforward application of the perturbation procedure would require 
very specific and lengthy calculations. 

It should also be noted that the modulations in the undular bore are not solely due to the 
dissipation. Rather, weak dissipation modifies the structure of the dissipationless undular 
bore. 



4 Modulation equations 

The derivation of the Whitham modulation equations for the KBB system (j^)) is based on 
the complete integrability of the unperturbed KB system 

ht + {hu)x + \u^xx = , . 
Ut + uux + hx = . 

That is, on the possibility to represent it as a compatibility condition of two linear equations 
for an auxiliary function ip: 

V-xx = Aij, = -\Bx^ + SVx (47) 

with ^ 

A={\~\^ -K S = - (a + , (48) 

where A is a spectral parameter. In the framework of this approach, the parameters A^ 
entering the periodic solution p5|) . ()16|) of Eqs. ()46|) have the following meaning. The 
second order differential equation (P7j) has two basis solutions and from which we 
can build the so-called 'squared basis function' 

g = i)^r- (49) 
It is easy to show that it satisfies the equation 

g^cxx - 2Axg - 4:Agx = 0, (50) 
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which after multipHcation hy g/2 can be integrated once to give 

\99x.-\9l-Ag^ = -P{\), (51) 

where the integration constant denoted by — P(A) can depend on the spectral parameter A. 
The second equation (jTTj) gives 

Qt = Bg^ - B^g. (52) 

In the finite-gap integration method (see, e.g. Kamchatnov, 2000), the periodic solutions 
are distinguished by the condition that -P(A) be a polynomial in A. Then g as a function of 
A should also be a polynomial in A. The one-phase periodic solution (fT5|l . (fT6|l corresponds 
to 

4 

P(A) = JJ(A - A,) = A^ - siX^ + S2\^ - S3A + S4 (53) 

i=l 

and 

g = \-^i. (54) 
Then from Eqs. ()5ip and ()52j) we find at once the relations ()15p as well as the equation for 

lie = 2^/P{ix) (55) 

whose integration yields (fT^ . 

As we see, the parameters Aj are the zeroes of the polynomial -P(A) which determine the 
periodic solution in the finite-gap integration method. At the same time, the parameters Aj 
are the most convenient modulation variables in terms of which the Whitham modulation 
equations assume the diagonal Riemann form (j2H) for the unperturbed KB equations (|46|1 
or its counterpart ()33|) for perturbed, KBB dynamics (j21). As was shown by Kamchatnov 
(2004), if the evolution equations are written symbolically as 

where the functions correspond to the "leading" , integrable part of the system, and the 
perturbation terms Rm can be slow functions of x and t and can also depend on the field 
variables u„ and their space derivatives, then the perturbed Whitham equations have the 
form 

where 

^. = -^^, . = 1,...,M. (58) 

Here the angle brackets denote averaging over the proper interval of x, M is the degree of 
the polynomial -P(A), is the order of the highest derivative Um""^ in A, and the index for 
the bracket means that A is put equal to Aj. 

In our case for the KBB system (0) we have M = 4 and 

N = 2: ui = h, U2 = u; 

R^ = 0, i?2 = i^Uxx, A2 = 0, dA/du = u/2- X. ^ ^ 
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Hence the perturbation terms on the right-hand side of Eqs. ()57p take the form 



For the one-phase modulated solutions of our present interest, one can replace averaging in 
(jFfjl . (jHHjl . (jHUjl with the averaging over the unperturbed periodic family Then, using 

diD), (jiH) we have 

(61) 



where the wavelength L is given by (j2!2j) . Further, taking account of the relations u = Si — 2/i, 
= 2^ P{p), Uxx = —AdP/djj, we have 

^\ 1 / f.-K-s,/2 =^#-1, (62) 



9/. 2L/ (A,-;x)v/pM " LdX, 



\ f / u\ dx 
{{X^ - u/2)guccx)i = J^J \^'~ 2/ ^^'"'''dil^'^ 

2 r {Xi-sj2 + fi){X,- fi)dP 



dfi 



(63) 



= ~j (V + \s,^^ - \s,X, + X^) ^^^^ci/z 
= -^j{^^-s,/A)^P(i^)d^i. 
Then, the characteristic velocities ()58|1 take the form 

coinciding with the unperturbed case ()2H) . while the perturbation terms in the Whitham 
equations (jHHjl are given by 



{dL/dx.)Y[^JX,-x~ ' 



A = TTTFTTTVTW^^ TT I \p - s,/ A)^!^) dii. (65) 



The integral here can be evaluated in terms of complete elliptic integrals. However, the 
resulting expression is very complicated and it is easier to deal with its unevaluated form 
(jnS)). The Whitham equations (jH^ . determine the evolution of the parameters Aj due to 
nonuniform modulation of the wave, and the small of effect viscosity. It is natural to expect 
that for the boundary conditions © the modulated wave will asymptotically, as t oo, 
tend to the steady undular bore solution described in Section 3. In the next section we shall 
find the corresponding stationary solution of the Whitham equations. 
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5 Steady solution of the Whitham equations 



We look for the solution of the Whitham equations ()33|) in the form 

Ai = Ai(^), e = x-ct, (66) 

so that we must have 

Motivated by the meaning of si/2 as the phase velocity one can suggest that equations (|F)7jl 
can be split in the following way: 



and 



where the factor 



^ = const (68) 



^ = 9. (69) 

"As 



Q = _^ _ s,/A)y^)d^^ (70) 

is the same for all z = 1, 2, 3, 4. 

For Eqs. (jUHj) . (jU^ to be consistent with Eq. Si must be an integral of equations 
(jnni)- In fact, we will show that the special structure of these equations provides actually 
three integrals Si, S2, S3. This statement can be proved with the use of the Jacobi identities 
(Jacobi 1884), which follow at once in the most convenient for us from the obvious identity 

trn,,A-Aj) ' ' 

where in the left-hand side we have a polynomial in A of the degree n — 1 which is equal to 
unity at n points A = Aj, i = 1, . . . , n, and hence is equal to unity identically. Then equating 
the coefficients of A™ at both sides of ()71|) we get n — 1 identities for m 7^ 0, 

V - = T = T ^'^ ' -0 (72) 

where prime means that all terms with the factor Ai are omitted in the corresponding sum, 
and the last identity for m = can be written in the form 



\n~l 



V 1 i-^y (73) 



where s„ = Yli In our case n = 4 and Eqs. and (f?^ yield 

^ = nV - = ^ = oy" ^ 



n,y.(A. - A,) ' de - ^ n,y.(A. - A,) 



d0 "^trn,y.(A.-A,) 



(74) 
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that is the system (jHH) has n — 1 = 3 integrals of motion 

Si = const, S2 = const, S3 = const. (75) 
Thus, in the steady solution only the last coefficient S4 varies with 6 according to the equation 
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where we have used the identity (ffH|) . 

Now, the zeroes A^, i = 1,2, 3, 4, are the solutions of the algebraic equation 

4 

P(A) = JJ(A - A,) = A^ - siA^ + S2X^ - S3A + S4 = (77) 

1=1 

ordered according to 

Ai < A2 < A3 < A4 (78) 

and for given si, S2, S3 they can be considered as known functions of S4. As a result, we 
arrive at the single first order differential equation 

dsA 8u /•^3(s4) 



de L 



/ {i^-Si/A)./P{ii}di^, (79) 



where L is given by Eq. For the solution under study, Aj = Aj(s4), hence in the right- 

hand side of (f7^ we have a known function of S4. The constants Si, S2, S3 can be expressed 
in terms of the initial parameters hi, h2, Ui, as in (jHHj) . To this end, we compare the equation 



ul = u^ - 4cM^ + 4(c^ - B)u^ + 8{Bc - A)u + const, (80) 
following from Eq. ()44|1 with i/ = 0, with the equation (see ()55|) ) 

fll = A{fi'^ - Sif/ + S2fJ? - S3/i + S4) . (81) 

Then taking account of the relation u = si — 2/i, these must coincide with each other, and 
so we find that 

si = 2c, S2 = c^-B, S3 = -{A + Bc), (82) 
where c is given by (|HH|) and (see 



To determine the interval within which the variable S4 can vary, we notice that at the 
leading and trailing edges of the undular bore the polynomial -P(A) has double roots, that 
is, its discriminant D vanishes. Hence, the limiting values of S4 must be the roots of the 
equation (see, e.g. Fricke 1924) 

D = gl- 21 gl = 0, (84) 
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where 

C/2 = S4 - IsiSs + j^sl, ^g^^ 



1 1 131212 

93 = QS2S4 + 4gSiS2S3 - 2J6S2 - YgSg - JqS^S4 

are invariants of the polynomial -P(A). Equation ()84|1 is cubic with respect to S4 and has 
three roots 54^^ < 54^'* < 54^'*. The variable n (and, hence, u) oscillates with finite amplitude 
as long as -P(A) has three real roots. Hence, S4 can vary between the two smaller zeroes of 
the discriminant D, 

sj^ < S4 < sf\ (86) 

Thus, all parameters in Eq. ()79p are completely determined and can be expressed in terms 
of hi, h2, Ui, so that dependence of S4 on 9 can be found by integration of Eq. (f7^ in the 
interval Eq. ()86|1 with the initial condition 

^ = s?^ at 6 = 00, (87) 

where we assume that the leading edge of the bore is located at 6* = 6*9. 
At the trailing edge Eq. (f?^ reduces approximately to 

^ = const ■(A3(54)-A2(^^4))'. (88) 

(1) 



Since in the vicinity of s\ we have 



then here 



which gives at once 



A2(s4), A3(s4) oc a/s4 - (89) 

"^^Cis,- .«) (90) 
S4 - s^4^ OC exp(C^) (91) 



where C is some constant proportional to u. Thus, we see that the trailing edge is formally 
located at ^ — oo, but with exponential accuracy we can take the width of the bore as 

const 

A ^ . (92) 

1/ 

The asymptotic analogous to Eq. has been obtained in (Gurevich and Pitaevskii 1987) 
and (Myint and Grimshaw 1995) for the KdV-Burgers equation. 

At the leading edge we have a soliton solution (fT^ . (fT^ with Aj = Xi{sf^). Its centre 
corresponds to /i(0) = A2 and hence Eqs. (fT^ give the values of the velocity Us and elevation 
hg at the centre of the leading soliton: 

Us = si- 2A2(sf ), K = Isj -S2- 2(A2(sf ))' + siA2(sf ). (93) 

The dependence of all values entering the right-hand parts of Eqs. (I93|) on the initial param- 
eters is given by Eqs. 
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Figure 8: The Riemann invariants as functions of S4. 

S4 




Figure 9: Dependence of S4 on 6 in the bore. 

To illustrate the developed theory, let us make some calculations and draw corresponding 
plots for specific parameters of the undular bore. Although the asymptotic approximation 
used in the derivation of the KBB system is consistent with the shallow water dynamics 
only for small values of the initial step /12 — /^i -C hi (see Eqs. (IHU jl - ipn]) ) it is instructive 
to consider the problem with noticeably distinct initial parameters hi and /12 for a better 
exposure of the details of the oscillatory structure. We choose u = 0.1 and 

ui = 0, hi = 1, h2 = 4. (94) 

Then we get 

U2 = 1.90, c = 2.53, A =2.53, 5 = 1.0 (95) 

and 

si = 5.06, S2 = 5.4, S3 = -5.06. (96) 
Equation (|H^ gives the limits for S4, 

= -7.74, = -6.4. (97) 

Solving Eq. (|77j) for Aj, we find the Riemann variables Aj as functions of S4 and the corre- 
sponding plot is shown in Fig. |H1 Integration of (ffn|) leads to dependence of S4 on 6 shown in 
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Figure 10: The Riemann invariants as functions of 6. 



Fig.ini Substitution of this dependence into Xii^s^), i = 1,2, 3, 4, yields the Riemann variables 
as functions of 9 depicted in Fig. [TUl As we see, the traihng edge is located at ^ — > — oo 
where X2{0) and A3 (6*) tend to the same limit A2(s4^'') = A3 (54^^). Finally, substitution of 
the slowly varying Riemann variables into (|15|) yields the profiles of velocity u{6) and water 
elevation h{6) in the bore; see Figs. ^2 and ^1 respectively. Obviously, the "camel hump" 
form of the lead soliton of the elevation profile in Fig. is due to the properties of the 
Kaup-Boussinesq system rather than actual properties of shallow water solitary waves. One 
should note, though, that within the range of applicability of the Boussinesq approxima- 
tion (i.e. for small enough initial steps) this deviation of the soliton shape from the regular 
shallow water soliton profile ceases to be visible. 



u(e) 




e 
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Figure 11: Velocity profile in the steady undular bore. 

One can see, that despite the different quantitative description, frictional and conservative 
shallow water undular bores are structurally similar in many respects (cf. El, Grimshaw & 
Pavlov 2001). There are, however, substantial qualitative differences. In particular: (i) the 
conservative undular bore expands in time while the frictional undular bore with the Burgers- 
like dissipative term asymptotically reaches a steady profile propagating with a single velocity 
c; (ii) the transition relations across the frictional and the dissipationless undular bores are 
different: for the frictional undular bore the transition relation (j38|) coincides with the jump 
condition following from the conservation laws ()34|). while in the dissipationless case (see 
(inH)) it coincides with the simple wave relation for the ideal shallow- water equations. 
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h(e) 




Figure 12: Elevation profile in the steady undular bore. 

We note in conclusion, that it is clear that the method of solution of the perturbed inte- 
grable Whitham equations used in this paper is essentially based on the special structure of 
the perturbation term (jU^ where the integral is actually a function of just a single parameter 
S4. Similar structure, however, can appear due to other than Burgers-type perturbation of 
the original integrable equation. Generally, this is the case for any form of perturbation term 
leading, after averaging, to the integrand in containing only symmetric functions of the 
Riemann variables Aj rather than individual Aj's. Then one can find sufficient number of 
integrals of the stationary Whitham equations to reduce the system to a single equation. 

6 Conclusions 

The formation of a shallow-water frictional undular bore has been studied analytically us- 
ing the Kaup-Boussinesq system modified by a small friction term. The main stages of the 
undular bore formation from the step-like initial profile were considered and the analytic so- 
lutions were constructed for the initial unsteady (dissipationless) and final steady (frictional) 
stages of the undular bore development, using the Whitham method. The perturbed inte- 
grable Whitham equations for the Kaup-Boussinesq-Burgers system were derived using the 
methods of finite-gap integration. It was shown that the stationary solution of the Whitham 
equations describing modulations in the steady undular bore is consistent with the jump 
conditions following from the original conservation laws for the KBB system. 

The theory developed in this paper shows that the Whitham method provides a general 
effective approach to describe frictional undular bores in perturbed integrable systems, and 
can be used in different physical contexts provided the dissipation is small enough to not 
prevent the generation of nonlinear dispersive waves, but sufficient to balance the combined 
action of nonlinearity and dispersion at large times. 
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